Experimental proposal for accurate determination of the phase relaxation time and 
testing a formation of thermalized non-equilibrated matter in highly excited quantum 

many-body systems 
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We estimate how accurate the phase relaxation time of quantum many-body systems can be 
determined from data on forward peaking of evaporating protons from a compound nucleus. The 
angular range and accuracy of the data needed for a reliable determination of the phase relaxation 
time are evaluated. The general method is applied to analyze the inelastic scattering of 18 MeV 
protons from Pt for which previously measured double differential cross sections for two angles in 
the evaporating domain of the spectra show a strong forward peaking. A new experiment for an 
improved determination of the phase relaxation time is proposed. The experiment is also highly 
desirable for an accurate test of a formation of thermalized non-equilibrated matter in quantum 
many-body systems. 

PACS numbers: 24.10.Cn; 24.60.Dr; 25.40.Ep 
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For some years one of us has emphasized that phase 
relaxation in a many-body system can be considerably 
longer than energy relaxation among independent parti- 
cle states U Q, 13 . Maybe the easiest experimental ac- 
cess to the problem can be found in low energy (10-80 
MeV) nucleon-nucleus scattering processes, where some 
of the oldest data have been available for half a cen- 
tury 0. It turns out p| that the phase relaxation time, 
Tnh = h//3 with P being the phase relaxation width 
HI yd - is considerably longer than the energy relaxation 
time, T erg = fi./r spr , obtained from standard estimates of 
the spreading width T spT of independent particle states. 

The significance of this problem results from the fact 
that an existence of long-living phase relations is of fun- 
damental importance in the study of relaxation phenom- 
ena in nuclear, atomic, molecular and mesoscopic many- 
body systems, and for many-qubit quantum computa- 
tion. In particular, if a phase relaxation time, which char- 
acterizes the lifetime of the "phase memory", is longer 
than the energy relaxation time, this effect could ex- 
tend the time for quantum computing d, [|| beyond the 
quantum-chaos border 0. 

The experimental result Q| showing that the intensity 
of the compound nucleus evaporation for = 60° ex- 
ceeds that for 8 — 150° by a factor of ~6 indicates that 
(3 ~ T T = H/t\ where is the compound nucleus decay 
width and is the compound nucleus life-time. How- 
ever, since the evaporation spectra were measured for 
two angles only, these data allow only a rough estimate, 
/3/r sp r ~ 10~ 5 Q- The purpose of this note is to de- 
termine whether more detailed experiments of the same 
type could improve this estimate. 

Our analysis for the improvement of the experimental 
determination of the phase relaxation time from nuclear 
evaporation data will be quite general and can readily 
be applied to any low energy nucleon-nucleus scattering 
showing forward peaking in the evaporation domain of 



the spectra. Yet, to be specific, we shall show numbers 
and graphs for inelastic scattering of 18 MeV protons 
from Pt, i.e. for improvement of one of the oldest exper- 
iments |4j- The questions we shall answer are: What is 
the best realistic accuracy of determination of the phase 
relaxation time from the data on forward peaking of evap- 
orating protons from the compound nucleus? What is the 
angular range and accuracy of data needed for a reliable 
determination of T p h? An advantage of the analysis is 
that it relies on relative values of the double differential 
cross sections which are usually determined experimen- 
tally with better accuracy than absolute cross sections. 

The proposed experiment is also desirable for an accu- 
rate test of the formation of thermalized non-equilibrated 
matter in compound processes. A manifestation of such 
a new form of matter, introduced in Refs. 13,11], would be 
equal slopes, i.e. nuclear "temperatures" g|, of the prop- 
erly scaled |4| proton evaporation spectra for forward and 
backward angles. Due to the insufficient statistics the 
data 4] indicate only approximate equality of the slopes 
for the forward and backward angles with about 20 per- 
cent uncertainty. We again point out that the proposed 
test only requires relative values of the proton emission 
intensities. 

The evolution of a nuclear reaction is usually consid- 
ered to proceed via a series of two-body nucleon-nucleon 
collisions, which successively form states of increasing 
complexity. On each stage of the reaction a distinction is 
made between continuum states and quasibound states. 
Emissions from the continuum states is related to mul- 
tistep direct reactions 0, IToL ITlT ] , and decay of the qua- 
sibound states originates multistep compound processes 
0,0]. The compound nucleus is formed at the last stage 
corresponding to the most complex configuration of the 
chain of quasibound states. The multistep direct reac- 
tions originate from the decay of the simplest configu- 
rations of the chain resulting in forward-peaked angular 
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distributions. In contrast, the multistep pre-compound 
and compound reactions are conventionally assumed to 
give rise to angular distributions symmetric about 90 de- 
grees. 

We use the exciton model to evaluate the rel- 
ative contributions of multistep direct, multistep pre- 
compound and compound nucleus processes for the p+Pt 
(E p = 18 MeV) inelastic scattering for the proton out- 
going energy of 7 and 6 MeV. Fitting the entire en- 
ergy range for forward angles we found |14| that, for 
the proton outgoing energy of 7 MeV, the compound nu- 
cleus cross section constitutes 90%, while multistep di- 
rect and multistep pre-compound are about 5% each. 
For the proton outgoing energy of 6 MeV, the compound 
nucleus cross section is about 98%. Therefore we ob- 
serve that even though the low energy 6-7 MeV outgoing 
proton spectrum is overwhelmingly dominated by com- 
pound reactions, the angular distribution is strongly for- 
ward peaked. Clearly a description of the decay of such 
thcrmalized but non-equilibrated matter requires a ma- 
jor modification of conventional theory of compound nu- 
cleus (see e.g. Kef. K. • originally formulated by Bohr, 
Bethe, Weisskopf, Wigner, Dyson and others. The basic 
assumption of the conventional theory is that thermal- 
ization of the compound nucleus guarantees a complete 
loss of memory of initial phase relations. A modification 
of this conventional picture of the compound nucleus was 
proposed by one of us in Refs. 0,0. The key element 
in the description of asymmetry of angular distributions 
around 90° cm. for evaporating particles is total spin 
off-diagonal correlation between compound nucleus par- 
tial width amplitudes. Such a correlation is neglected in 
a conventional picture of compound nucleus. Following 
0,13 we have 
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where overlines denote ensemble averaging. Here J\ J 2 
are the compound nucleus total spin values, 7Ti,7T2 are 
parity values, E^ are resonance energies with n being 
running indices, and D is average level spacing of the 
compound nucleus. The a(b) indices specify the orbital 
momenta l ai2 {h li2 ) , the channel spins j ai 2 (i6 12 ), and 
the microstates a(b) of the target nucleus and residual 
nucleus, respectively. Accordingly, a\ — a 2 denote the 
ground state of the target, and b\ — b 2 specify the mi- 
crostates of the residual nucleus. The above correlation 
between the partial width amplitudes leads to a corre- 
lation between fluctuating compound nucleus ^-matrix 
elements carrying different total spin values: 
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Here S'l b (E) are compound nucleus S- matrix elements 
with total spin J and the brackets (...) denote the en- 
ergy E averaging. For finite values of j3/T\ nonvanish- 
ing of the spin off-diagonal correlations in Eq. (2) reflects 
nonvanishing of the interference between resonance levels 
with different total spins upon the energy averaging. 

For the correlation between S'-matrix elements carry- 
ing the same total spin and parity values and the same 
microstates a± = a 2 and b\ = 62 but different orbital 
momenta and/or channel spins we have 

(s j a : bl (Eys j a : b2 (E)) = [i ( E ) 1 2 } (fe 2 i E ) i 2 )] 1/2 - 

(3) 

The above equation results from a strong correlation be- 
tween partial width amplitudes -y^ 01 ^ 1 ) anc ] ^■^ a ^ b2 > 
with ai = a 2 and h = b 2 but l ai ^ l a2 , l bl ^ l b . 2 , 
j ai io 2 , jb t 7^ jb 2 ■ Such a correlation is referred 
to 0, Q as the continuum correlation. Note that such 
a strong correlation between reduced width amplitudes 
corresponding to the same total spin and parity values 
but different orbital momenta was experimentally re- 
vealed for a number of compound nuclei in the regime 
of isolated resonances . 

For f3 ^> T\ the spin off-diagonal correlations in Eq. 
(2) result in the angular distributions symmetric around 
90° cm. recovering a conventional picture of compound 
nucleus. However, if (3 < , i.e the phase relaxation 
time T p h is comparable or longer than the average life- 
time of the compound nucleus, this allows us to describe 
a strong asymmetry of the angular distributions around 
90° cm. of the evaporating yield. 

For the treatment presented here we follow Ref. 
We neglect the intrinsic spins of the scattering partners in 
the entrance channel and proton spin in the exit channel. 
Since evaporated protons carrying orbital momenta I > 1 
are significantly sub-barrier due to both the centrifugal 
and Coulomb barriers, we take TJ>i = 0, where XJ are the 
transmission coefficients for the inverse process of capture 
of the proton by the residual nucleus. Then the double- 
differential cross section in the evaporation domain of the 
spectra has the form 
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Here, a(e) is the angle-integrated cross section for the 
evaporation of a proton with the energy e, and Pl{cos9) 
denote the Legendre polynomials of order L. The coeffi- 
cients of the angular decomposition are given by 
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FIG. 1: Dependence of Ti/T on /3/r T obtained from Eq. 
(6) for IZi = 6 (solid line), 5 (dashed line), and 7 (dotted 
line). The shadowed area shows the range of possible Ti/To 
and /3/r T which are in accordance with the experimental data 
of Ref. |4J . The dots marked with numbers on the solid line 

correspond to 1: (&■, ^)=(0, 0.72); 2: (0.1, 0.63); 3: (0.2, 

0.52); 4: (0.3, 0.39); 5: (0.367, 0.231); 6: (0.3, 0.13); 7: (0, 
0.05). 



One can see that if this phase memory time r p h = H//3 is 
about or longer than the average life-time t* = h/T^ of 
the compound nucleus the evaporation yield is emitted 
asymmetrically about 90° cm., i.e. the memory about 
the direction of the initial beam remains. However, if the 
phase memory time is much shorter than the average life- 
time of the compound nucleus then the spin off-diagonal 
correlations vanish, the memory on the direction of the 
initial beam is lost, and an isotropic angular distribution 
around 90° cm. is obtained. 

We now turn to the experimental data available from 
Ref. in order to estimate bounds for the possible val- 
ues of p. Since in Ref. only relative yields for two 
scattering angles are reported we focus on the analysis of 
the shape of the angular distributions 1(6) = KcPa/dflds 
without paying attention to the angle and the energy in- 
dependent prefactor K. 

For e — 7 MeV, we find from the experimental data a 
ratio 1(6 = 60°)//(6> = 150°) « 6. Thus, we consider the 
relation 
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where 1(6) is given by Eqs. (4) and (5) and depends on 
the two parameters P/T* and Ti/T . Therefore, Eq. © 
implicitly defines a curve in the parameter space of P/T* 
and Ti/To, which is plotted in Fig. ^ We can either 
solve Eq. © for p/T* as a function of Ti/To or vice 
versa. In any case we have an underdetermined equa- 
tion demonstrating that measurements of the evapora- 
tion yields for two angles only does not allow an un- 
ambiguous determination of the phase relaxation time. 
Indeed, in Fig. ^ any of the seven dots as well as any 
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FIG. 2: Angular distributions for the different sets of P/T* 
and Ti/To as specified in Fig. 1. In the inset a magnification 
is shown for the angular range we used for our analysis (see 
text). 



point on the solid line corresponds to IZi = 6. Yet, mov- 
ing along the solid line, P/T' changes from to 0.367. 
Since the value of Ti/T can be obtained from model 
calculations this clearly demonstrates that measurement 
of the evaporation yields for two angles only permits to 
accurately determine only the upper experimental limit, 
P/T* — 0.367. Estimating the total decay width from 
the systematics in Fig. 7 of Ref. ^3] we obtain T* = 0.02 
keV and, therefore,/? < 7 eV. Taking into account that 
T spr ~ 1.5 MeV |5[ we observe that the phase relax- 
ation time is at least five orders of magnitude longer than 
the energy relaxation time. Yet, it is still about fifteen 
orders of magnitude shorter than the Heisenberg time, 
h/D, where the average level spacing of the compound 
nucleus, D ~ 10~ 20 MeV, has been estimated using the 
Fermi-gas model [T^ |. 

Due to the uncertainties of the data reported in Ref. 
[5J, we also included curves of Ti/T and P/T" in Fig. [I] 
resulting from Eq. © but with its right hand side being 
5 and 7. 

For a more accurate determination of the p, rather 
than just an estimation of its upper limit, one has to 
measure the evaporation yields for more than two angles. 
Therefore, we analyze the sensitivity of the shapes of the 
angular distributions for different sets of P/T' and T\/Tq 
belonging to the manifold obtained from the solution of 
Eq. (6). 

In Fig. 2 we present the angular distributions for the 
seven sets of P/T* and Ti/T marked with dots and cor- 
responding numbers in Fig. 1. The angular distribu- 
tions are normalized in such a way that 1(6 = 60°) = 6 
and 1(6 = 150°) = 1 for each of these curves. One 
does observe that the angular distributions change ap- 
preciably with the change of P/T' and T\/T§ values. To 
quantify this sensitivity we determine 6 op t for which the 
ratio I (6 pt) /I (6 = 170°) is most sensitive to different 
P/T* and Ti/To values. We find that, for Kj = 6,5,7, 
6 opt = 118°, 116.3°, 119.4°, respectively. 

The dependence of I(6 opt )/I(6 = 170°) on P/T* is pre- 
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FIG. 3: The dependencies of I(6 opt )/I(6 = 170°) on /?/T T . 
The solid line is obtained for IZi = 6 (see Eq. JJJ) the dashed 
line for IZi = 5, and dotted line IZi — 7. Dots marked with 
numbers on the the solid line correspond to the /3/F T and 
Ti/To values specified in Fig. 1. 



sented in Fig. 3 with the solid line corresponding to the 
ratio TZi = 6. We find that for values of I(9 op t)/I(8 = 
170°) < 2.7 - 2.8 the ratio I(6 opt )/I{9 = 170°) is rather 
sensitive to the 0/T\ as can be seen in the lower part 
of Fig. |21 If the experimental value lies in this range, 
an accuracy of I(9 opt )/I(0 — 170°) of about 5% al- 
lows a determination of (3/T^ with a minimal uncer- 
tainty of about 10%. For a too low experimental value 
of I(0 op t)/I(9 = 170°), say less than 0.84, one can deter- 
mine only the upper limit of (3/T^ < 0.02. 

On the other hand, in the upper part of the curve, for 
I {0 opt) 1 1 {0 = 170°) > 2.7 - 2.8, the ratio I(0 opt )/I(6 = 
170°) shows a rather weak dependence on 0/T'. In par- 
ticular, for an assumed accuracy of the data of about 
5%, the value of f3/T^ can be determined with an uncer- 
tainty of only ~ 50 — 100%. If the experimental value of 
I (6 opt) I '1(0 = 170°) is larger than 2.85 one can determine 
only the upper limit, (3/T^ < 0.2. 

Similar conclusions about the sensitivity of 
I [0 opt )/ 1(0 = 170°) to /3/r T can be drawn from 
Fig. 3 for IZi = 5 and 7. According to Fig. [3] an 
improvement of the accuracy might be found if the data 
for forward scattering (9 < 30°) would also be available. 

It should be noted that a manifestation of a formation 
of the thermalized non-equilibrated matter has been also 
identified from a strong asymmetry around 90° cm. of 
evaporating protons in the Bi(7,p) photonuclear reaction, 
see [y and references therein. Other examples of a strong 
forward peaking of evaporating protons are found, e.g., 
in Bi(p,p') and Bi(n,p) processes with the 62 MeV en- 
ergy of initial beam, see |5j and references therein. Even 
though the later examples do demonstrate a formation 
of the thermalized non-equilibrated matter in compound 
processes, a determination of the phase relaxation time 
from these data is not unambigous since, for such high 
energy of the initial beam, there is a high probablility for 



a second and third chance proton evaporation. 

The possibility that in highly excited many-body sys- 
tems the phase relaxation can be much longer than the 
energy relaxation may have significant implications for 
quantum computing 0, |(| as well as, e.g., time-delayed 
"statistical" ionization of many-electron quantum dots 
and atomic clusters (see, e.g., [jjj and references therein). 
A possible presence of the effect of anomalously slow 
phase relaxation [2(| in chemical reactions (see [2(j and 
references therein) would require a modification of the 
statistical theories - phase space and transition state the- 
ories (see, e.g., [2l| and references therein). Yet, the nu- 
clear data indicating an existence of anomalously slow 
phase relaxation, which is much slower than the energy 
relaxation, are largely unrecognized by nuclear physicists 
and unknown outside the nuclear physics community. 
In many fields, including statistical physics, the notion 
"thermalization" or "energy equilibration" is considered 
to be equivalent to the notion "statistical equilibrium" . 
This note is a step towards changing this undesirable sit- 
uation. 

In conclusion we have proposed a general method to 
estimate the accuracy of the determination of the phase 
relaxation time from data on forward peaking of evap- 
orating protons from compound nucleus. The angular 
range and accuracy of the data needed for a reliable de- 
termination of the phase relaxation time have been eval- 
uated. The general method has been applied to the anal- 
ysis of inelastic scattering of 18 MeV protons from Pt 
for which previously measured double differential cross 
sections for two angles in the evaporating domain of the 
spectra show a strong forward peaking. We found that 
a new measurement of the angular distributions of evap- 
orating protons in the Pt(p,p') inelastic scattering for a 
wider angular range should permit an accurate determi- 
nation of the phase relaxation time. The experiment is 
also highly desirable for an accurate test of the formation 
of thermalized non-equilibrated matter in compound pro- 
cesses. Our analysis for the improvement of the experi- 
mental determination of the phase relaxation time from 
nuclear evaporation data can readily be applied to any 
low energy nucleon-nucleus scattering showing forward 
peaking in the evaporation domain of the spectra. 
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